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■^ I Abstract 

!>■. 

O , We discuss the phase structure of A^ D4 branes wrapped on a tempo- 

ral (Euclidean) and a spatial circle, in terms of the near-horizon geometries. 
This system has been studied previously to understand four dimensional pure 
/S ' SU{N) Yang-Mills theory (YM4) through holography. In the usual treat- 

5r ■ ment of the subject, the phase transition between the solitonic D4 brane 

and the black D4 brane is interpreted as the strong coupling continuation 
of the confinement/deconfinement transition in YM4. We show that this 
interpretation is not valid, since the black D4 brane and the deconfinement 
phase of YM4 have different realizations of the Zn centre symmetry and 
cannot be identified. We propose an alternative gravity dual of the confine- 
ment/deconfinement transition in terms of a Gregory-Lafiamme transition of 
the soliton in the IIB frame, where the strong coupling continuation of the 
deconfinement phase of YM4 is a localized D3 soliton. Our proposal offers 
a new explanation of several aspects of the thermodynamics of holographic 
QCD. As an example, we show a new mechanism of chiral symmetry restora- 
tion in the Sakai-Sugimoto model. The issues discussed in this paper pertain 
to gravity duals of non-supersymmetric gauge theories in general. 
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1 Introduction 

Strongly coupled field theories such as QCD and various condensed matter 
systems are not amenable to perturbative calculations and require numer- 
ical or other non-perturbative tools. A powerful new tool in this context 
is holography, which maps gauge theories to gravity in a higher dimension 
[H El El in El E] • Holography, however, is most reliable and quantitative for 
special supersymmetric gauge theories. In the generic case, including for pure 
Yang Mills theories, inferences from gravity tend to be less quantitative due 
to various degrees of extrapolations, although such exercises have proved to 
be rather valuable, e.g., in the context of holographic QCD O [6l [71 [H [9] . 

Recently, it was pointed out in [10] that there are problems with a naive 
application of holography to a two dimensional large A^ bosonic gauge the- 
ory at finite temperature. The dual of this system is based on gravity back- 
grounds involving A^ D2 branes on T^. The phase diagram of gravity, inter- 
preted naively, does not admit a unique continuation to the regime of gauge 
theory, owing to its dependence on the boundary condition for fermions on 
the brane. However, it was shown in [10] that the phase structures in grav- 
ity and gauge theory could be smoothly connected through an appropriate 
choice of fermion boundary conditions. 

In this paper, we will show that similar problems arise when we holo- 
graphically analyse more general non-supersymmetric gauge theories at finite 
temperature. The particular example we will focus on is that of holographic 
QCD from D4 branes [SI El [71 El [9)^1 . We will discuss some problems with 
the usual correspondence between the confinement /deconfinement transition 
in QCD and the Scherk-Schwarz transition between a solitonic D4 brane 
and a black D4 brane in the gravity dual. Some of these problems were 
first discussed in [TT] (see section [3?2] for further details). Especially we will 
show that the black D4 brane cannot be identified with the (strong coupling 



^We will briefly discuss the generalization to a p dimensional bosonic gauge theory and 
its holographic dual in section [T] 



continuation of the) deconfinement phase in QCD in four dimensions. As 
a resolution of these problems, we will propose an alternative scenario in 
which the confinement /deconfinement transition corresponds to a Gregory- 
Laflamme transition [121 [131 [S] between a uniformly distributed D3 soliton 
and a localized D3 soliton in the IIB frame. 

The scenario we propose suggests that we need to reconsider several pre- 
vious results in holographic QCD including the Sakai-Sugimoto model [9]. 
One important ingredient in the Sakai-Sugimoto model is the mechanism of 
chiral symmetry restoration at high temperatures [11]. We will propose a 
new mechanism for chiral symmetry restoration in our framework. 

The plan of the paper is as follows. 

Section |2] contains a short review of finite temperature QCD and the 
holographic approach to four dimensional Yang Mills theory (YM4)o using D4 
branes. In Section 13.11 we discuss the gravity theory at a finite temperature 
and a 'Scherk-Schwarz' transition between a solitonic D4 brane and a black 
D4 brane. In Section 13.21 we recall the conventional correspondence between 
this transition and the confinement /deconfinement transition of the Yang 
Mills theory, and discuss problems with this correspondence, with special 
emphasis on the mismatch of Zn x Zn centre symmetry between the black 



D4 phase and the deconfined 4 dimensional YM theory. In Section 14. H 
we discuss the phase structure of the D4 branes using an unconventional 
(periodic) fermion boundary condition along the Euclidean time circle and 
show that the high temperature in gravity is described by a localized solitonic 
D3 brane (or its T-dual) whose centre symmetry precisely matches with the 
deconfinement phase of four dimensional YM theory. This leads us, in Section 
14.21 to propose a new representation of the confinement /deconfinement phase 
transition in terms of a GL transition between the solitonic D4 and the (T- 
dual of) the localized solitonic D3 (see Table |2|). In Section [5l we show some 



^ Throughout this paper we will often use the notation YM4 for four dimensional Yang 
Mills theory and SYM5 for five dimensional super Yang Mills theory. 



new correspondences between phenomena in QCD and their counterparts in 
gravity following our proposal. In Section |6l we suggest a new mechanism 
of chiral symmetry restoration in the Sakai-Sugimoto model in keeping with 
our proposal. Section [7] contains the concluding remarks and some open 
problems. 

2 Review of QCD and its gravity dual 

2.1 Finite temperature QCD 

In this section, we briefly review some salient properties of four dimensional 
SU{N) pure Yang-Mills theory at a finite temperature at large A^, a system 
which we will subsequently investigate through holography. 

An important symmetry in this theory is the Zj^ symmetry along the 
temporal cycle, which is the centre of the SU{N) symmetry. The order 
parameter of this symmetry is the temporal Polyakov loop operator 

Wo = —TTPe'^o^odx°. (1) 

A^ V / 

If this operator is zero, the Z^v symmetry is preserved and, if not, it is 
broken. Especially, if the Zj^f symmetry is preserved, physical quantities 
do not depend on the temporal radius /3/27r at 0{N'^) order due to 'large 
A^ volume independence' [I5l [161 [E]- If the gauge theory is on a torus 
S\^ X S\^ X ■ ■ ■ , the large A^ volume independence is generalized and physical 
quantities do not depend on Lj, if the Zn symmetry along the i-th direction 
is preserved [T5] . 

Let us consider the phase structure of the Yang-Mills theory. At low tem- 
peratures, the confinement phase is the dominant thermodynamic phase. In 
this phase, the Z^ symmetry is preserved {Wq = 0). At high temperatures, 
the deconfinement phase, in which the Z^ symmetry is broken {Wq ^ 0), 
dominates. There is a confinement/deconfinement transition at an interme- 
diate temperature. From studies in large A^ lattice gauge theories, this is 



expected to be a first order phase transition [13, [201 ED 122] • 

As we remarked above, the free energy in the confinement phase does not 
depend on the temperature through the 'large TV volume independence'. It 
implies that the entropy becomes zero at 0{N'^). On the other hand, the 
entropy in the deconfinement phase would be 0{N'^). This has a simple 
physical interpretation: in the confinement phase, the spectrum excludes 
gluon states, which have 0{N'^) degree of the freedom, and consists only of 
gauge singlet states like glueballs, leading to an 0(1) entropy. 

2.2 Holographic QCD 

In this section, we will review the construction of four dimensional SU{N) 
pure Yang-Mills theory from A^ D4 branes [S]. Let us first consider a 10 di- 
mensional Euclidean spacetime with an S^ and consider D4 branes wrapping 
on the 5*^. We define the coordinate along this S^ as x^ and its periodicity 
as L4. The effective theory on this brane is a 5 dimensional supersymmetric 
Yang-Mills theory on the S}^^. For the fermions on the brane, the boundary 
condition along the circle can be AP (antiperiodic) or P (periodic); to specify 
the theory, we must pick one of these two boundary conditional. Let us take 
the AP boundary condition. This gives rise to fermion masses proportional 
to the Kaluza-Klein scale I/L4, leading to supersymmetry breaking (this 
is called the SS — Scherk-Schwarz — mechanism). This, in turn, induces 
masses for the adjoint scalars and for A4, which are proportional to Xa/L^ at 
one-loop. Therefore, if A4 is sufficiently small and the dynamical scale Aym 



^ Defining a spinor field on a manifold M requires specifying a "spin structure" , which 
means specifying a consistent choice of the sign of the spinorial transformation (recall that, 
given a Lorentz transformation A on vectors, the corresponding spinor transformation S'(A) 
is ambiguous up to a sign S'(A) o — ^(A)). On a simply connected Euclidean manifold M 
there is a unique spin structure; all cycles are contractible and the fermion b.c. along any 
such cycle is antiperiodic. If M has the topology S^ x N where N is simply connected 
(generalizations to multiple S^ factors are obvious), then the spin structure on M has a 
two-fold ambiguity, where a spinor can be either periodic (P) or anti-periodic (AP) along 
the S^. See, e.g. [13]. For examples relevant to this paper, see [S]). 



is much less than both the above mass scaleqj, then the fermions and adjoint 
scalars are decoupled and the 5 dimensional supersymmetric Yang- Mills the- 
ory is reduced to a four dimensional pure Yang-Mills theorja. 

By taking the large A^ limit of this system a la Maldacena at low temper- 
atures, we obtain the dual gravity description of the compactified 5 dimen- 
sional SYM theory [5], which consists, at low temperatures, of a solitonic D4 
brane solution wrapping the S^^. The explicit metric is given by 



ds =a' 



u 



VdiX^ 



dt' + V dx^ + Mu)dxl I + ^^ r^ + u'dni] 



There is also a non-trivial value of the five form potential which we do not 
write explicitly. Here A5 = A4L4 is the YM coupling on the D4 world-volume 
and di is found by putting p = 4 in the general formula 

dp = 27-2p^(9-3p)/2r n-^\ . (3) 

Since the X4-cycle shrinks to zero at m = mq, in order to avoid possible 
conical singularities we must choose the asymptotic radius L4 as follows 

-L4 yd^X^ _i/2 ,.s 

^ = ^^—"0 • (4) 

With this choice, the contractible a;4-cycle, together with the radial direction 
M, forms a so-called 'cigar' geometry which is topologically a disc. In order 
that the fermions are well-defined on this geometry, they must obey the AP 



^ Recall that Kym = Aexp[— 1/(6A4)], where A4 is the value of the coupling at some 
cut-ofF scale A and h is the first coefficient of the /3-function; if A4 is sufficiently small it 
is possible to ensure Kym ^ A4/L4, 1/L4 since as A4 goes to zero, exp[— 1/(&A4)] goes to 
zero faster than A4. 

^In this construction of the 4 dimensional Yang-Mills theory from the 5 dimensional 
SYM through the Kaluza-Klein reduction, we have implicitly assumed that the long string 
modes are suppressed. We will come back to this problem later on. 



boundary condition. The AP boundary condition along X4 is, of course, 

consistent with the choice of fermion b.c. in the boundary theory. 

The leading order gravity solution described above is not always valid. 

E.g. in order that the stringy modes can be ignored, we should ensure that 

the curvature in string units must be small. In other words, the typical length 

/ , ,„\ 1/2 

scale of this solution near u = uo, viz, I = la'y/d^X^UQ j , must satisfy 

/ ^ Va' p]. This condition turns out to be equivalent to 

A4 > 1. (5) 

This is opposite to the condition which we found for the validity of the four- 
dimensional gauge theory description (see footnote H]). Thus, this gravity 
solution can describe the 5 dimensional SYM but cannot directly describe 
the 4 dimensional YM theory, inferences about which can only be obtained 
through extrapolation. This is a common problem in the construction of 
holographic duals of non-supersymmetric gauge theories. This has been dis- 
cussed at length; in particular, the gravity description, which necessitates 
extrapolation to strong coupling, has been likened (cf. [7] and [6], p. 196- 
197) to strong coupling lattice gauge theory. Many interesting results, in- 
cluding the qualitative predictions in [SI El [TJ [9] , have been obtained using 
this prescription. 

Let us explore some properties of the soliton solution ([2]). This solution is 
expected, on the basis of several arguments, to correspond to the confinement 
phase in the four dimensional gauge theory. For example, it can be shown 
that the classical action of this solution at a finite temperature 1//3 is [H 

ElEsfl 



Ssda/N^Vs = -C4A4/3LJ^ (6) 

where V3 is the volume of the three spatial dimensions (xi, X2, X3) and C4 is a 



^In Scction lOl we will describe more details of holographic QCD at finite temperatures. 
Eq. (jni) describes the stable phase of the system only at a sufficiently low temperature. 



constant which is obtainable from the formula below [21 [TU] by putting p = 4 



^^ 2ii-2p^(i3-3p)/2r((9-p)/2)af-^)/(^-^)' """ And'/' ^^^ 



Thus the free energy of this solution, F = Ssda/ f!^, is independent of tempera- 
ture, and hence the entropy is zero at 0{N'^) order. These facts are consistent 
with the interpretation of this solution as the confinement phase in the large 
A^ gauge theory, in which the temporal Z^ symmetry is preserved. 

Note that the action ([6]) divided by L4 does depend on L4, implying 
that the Zj^ symmetry along the X4-cycle is broken. Although this observa- 
tion might appear irrelevant for the 4 dimensional Yang-Mills theory, it is 
this symmetry breaking which allows a conventional KK reduction to four 
dimensions (see the next two sections for more details). 

3 Finite temperature holographic QCD 

3.1 Confinement /deconfinement transition in Holo- 
graphic QCD in the standard scenario 

In existing holographic studies of the thermodynamics in the Yang-Mills the- 
ory, it is assumed that the black D4 brane solution corresponds to the decon- 
finement phase and the Scherk-Schwarz transition in the gravity is related 
to the confinement/deconfinement transition. In this section, we review this 
scenario and we will explain its problems in section 13.21 

To discuss holographic QCD at finite temperatures, we begin by com- 
pactifying Euclidean time in the boundary theory on a circle with period- 
icity /3 = l/T. To permit description in terms of four dimensional gauge 
theory, the temperature should, of course, be kept far below the KK scale, 
T <^ A4/L4 y\. In order to determine the gravitational theory, we need to 



^The precise conditions for the decoupling of KK modes from YM4 are T <C 1/^4 for 
the fermions and T <C A4/L4 for the bosonic modes (with A4 ^ 1). Since the second 
condition includes the first one at weak coupling, we omit the first condition. 



fix the periodicity of fermions in tlie gauge tlieory along tlie time cycle. In 
the existing literature, the fermion boundary condition is chosen to be AP 
according to the usual practice for thermal fermions (we will see in section H] 
that this choice is not imperative here). The resulting fermionic bc'qj: (AP, 
AP) along (t, X4) lead to a Z2 symmetry of the system under t -H- X4. In the 
gravity dual, we must therefore, include the Z2-transform of the solution (|2]), 
which is a black D4 solution: 



ds =0L 



u 



\/d\\ 



5 



j,{u)dt^ + V dx\ + dxA + ^^ f ^ + uHVl\ 



no\3 



Mu) = 1 - ( ^j • (8) 



Here Uq in this metric is related to /3 as in 

-uo'''. (9) 



/3 y/ckM -1/2 



2tt 3 

This solution has a contractible t-cycle and the fermions must obey the anti- 
periodic boundary condition (see footnote [3]) , consistent with the AP b.c. in 
the gauge theory. Analogously to (l5j), the gravity description of this solution 
is reliable if A4 > fS/L^. 

The classical action density can be evaluated, as before, yielding 

Sbd4/N^Vs = -CX^Llr'- (10) 

In contrast with ([6]), the free energy now is a function of /3 and hence the 
entropy is 0{N'^) which is appropriate for a description of the deconfinement 
phase in the gauge theory. 

By comparing (|6]) and (fTOj) . we see that at low temperatures (large /3) the 
solitonic D4 solution dominates, while at high temperatures the black brane 



* Here and elsewhere the boundary conditions will always refer to those of the boundary 
theory along S^ x S\^, respectively. 

^In several equations for the black D4 brane in this paper, we use the four dimensional 
gauge coupling A4 = A5/L4, which obscures the /3 o L4 symmetry. If we write these 
equations using A5 instead of A4, we will see that they are rather simply obtained from 
the corresponding equations for the solitonic D4 brane via the exchange /? O L4. 

9 



dominates. The transition between the two solutions, which we will call the 
"Scherk-Schwarz (SS) transition", occurs at 



(3ss = U. (11) 

In the standard discussions of holographic QCD [21 [S], [TT], this transition is 
interpreted as the (continuation of the) confinement /deconfinement transi- 
tion in the gauge theory. In the following sections, we will explain various 
problems with this interpretation and present a resolution by proposing an 
alternative scenario. 

3.2 Problems with the standard gravity analysis 

The previous subsection describes the usual relation between the confine- 
ment /deconfinement transition in 4 dimensional Yang-Mills theory and the 
SS transition in holographic QCD. However as we will now discuss, there 
are several problema^i associated with this correspondence. Some of these 
problems were mentioned earlier in [TT] . 

In order to understand these problems, it is convenient to consider first 
the phase structure of the 5 dimensional SYM on the S\ x S\^. See Figure 
[U As we have seen, the gravity description is valid in the strong coupling 
regime of this theory (the blue region in Figure [T]) which is characterized by 
the solitonic D4 solution and the black D4 solution. On the other hand, the 
4 dimensional Yang-Mills theory is realized at weak coupling (A4 <^ 1) and 
low temperature (/3 ^ L^^/X^) (the upper green region in Figure [I]) which is 
characterized by the confinement phase and the deconfinement phase. The 
usual proposal is that the solitonic D4 and black D4 solutions correspond to 
the confinement and deconfinement phases respectively. 

To examine this proposal in detail, it is useful to consider the realization 
of the Ztv X Z]^ centre symmetry in various phases of the five-dimensional 



^"Thc problems discussed in this section are not about the use of holography in pre- 
dicting the phase diagram of strongly coupled 5 dimensional SYM, but rather about its 
extrapolation to four dimensional Yang-Mills theory. 

10 
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(A4 » 1) 
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(A4 > /3/L4) 

^^4 = 



I 



A4 



Figure 1: The phase structure of the five dimensional SYM (SYM5) on 5^ x 
S\^ with the (AP,AP) boundary condition. The gravity analysis is valid 
in the strong coupling region (the blue region). The 4 dimensional YM 
description is valid in the upper green regime. The lower green regime is the 
mirror of the upper one via the Z2 symmetry /3 <H- L4. The solid red line in 
the green regions depicts the known confinement/deconfinement transitions 
in gauge theory and its Z2 mirror. The solid red line in the blue region 
is the SS transition (fTTjl . These lines demarcate phases characterized by 
the Zn X Zn order parameters Wq and W^ (see Table [1]) . The solid black 
lines correspond to a minimal extrapolation of the phase boundaries through 
the intermediate region. The dotted line denotes another possible phase 
transition which is allowed by the Z2 symmetry but across which the Zat x Ztv 
symmetry does not change. 
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'Confinement' phase in YM4 

'Deconfinement' phase in YM4 

Z2 mirror of 'confinement' phase in YM4 

Z2 mirror of 'deconfinement' phase in YM4 

Sohtonic D4 phase 

Black D4 phase 

Locahzed sohtonic D3 phase 



1^0 = 0, 






7^0 


7^0 


= 


7^0 


7^0 


= 


7^0 



AP/P 

AP/P 

AP 

AP 

AP/P 

AP 

P 



Table 1: The values of the Polyakov loops in some phases of the 5 dimensional 
SYM (SYM5) on S*^ x S\^ (focussing on features related to the physics of the 
four dimensional YM theory). The boundary condition for the fermion along 
S\ is shown in the right column. (We have fixed the AP b.c. along the 5*^^.) 
The 'confinement' phase in YM4 here refers to the phase of SYM5 whose KK 
compactification, at small enough L4, coincides with the usual confinement 
phase of YM4. Similarly, the 'deconfinement' phase in YM4 here refers to the 
phase of SYM5 whose KK compactification, at small enough L4, coincides 
with the usual deconfinement phase of YM4. The phases labelled by gravity 
solutions, valid at suitable parameter regimes, denote the respective phases 
of SYM5 at those parameter regimes. 



SYM theory along t and X4. The Z^v symmetry along t is characterized by 
the order parameter Wq which is the temporal Polyakov loop ([T]). Similarly 
the Z]si symmetry along 0:4 is characterized by the order parameter W4 which 
is the Polyakov loop around the X4-cycle: 



1^4 = — TrPe^^o^'^*'^^' 

N 



(12) 



In Table [T] we have collected values of {Wq, W^ for the phases of 5D SYM 






discussed in this papeio. 

Before testing the above proposal for holographic QCD in terms of these 
order parameters, let us see how to obtain this table. The 'confinement' phase 



^"'^One cannot, a priori, rule out additional phases in the intermediate coupling regime 
(the white region in Figure [1] and [3]). However, the lattice super Yang Mills study in [53] 
has analyzed such intermediate coupling regimes in the context of Dl branes compactified 
on a periodic circle, and has not found any novel phases. This result may indicate that 
there are no other phases in the present case with D4 branes as well. 
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denotes a phase of the 5D SYM whose KK compactification, at small enough 
L4, coincides with the usual confinement phase of YM4. The latter, clearly 
has Wo = 0. Now, the fact that we have obtained the 4D phase through a KK 
reduction, along L4, implies that the corresponding phase of the 5D theory 
must have W4 ^ 0. This follows from the fact [T71 |2Z] that in a W4 = 
phase, KK reduction does not work since the effective KK scale in a large A^ 
gauge theory becomes l/i^NL^j in stead of the usual I/L4, and goes to zero 
at large iVo. Thus, the 'confinement' phase must have Wq = 0, W4 7^ 0, as 
shown in the first row of Table [1] By similar arguments, the 'deconfinement' 
phase of 5D SYM which coincides with the deconfinement phase of YM4 
must have W4 7^ (for validity of the KK reduction) and Wq 7^ (to exhibit 
4D deconfinement). This is shown in the second row of Tabled] 

The next two rows in Table [1] describe two additional phases of SYM5, 
which appear only for the (AP,AP) b.c. via the Z2 symmetry. The third row 
describes the Z2 mirror of the confinement phase of YM4, which is character- 
ized by Wq 7^ and W4 = 0. Since W4 = 0, this phase is a fully 5 dimensional 
object and is not related to the YM4 (indeed, so far as the Polyakov loop 
order parameters are concerned, this phase can be identified with a decon- 
finement phase of SYM5). Another phase is the mirror of the deconfinement 
phase in YM4, which is characterized by Wq 7^ and W4 7^ 0. Since this 
phase has the same order parameters as the original deconfinement phase, it 
is not obvious whether these two phases are distinct or smoothly connected 
through a cross over (we represent the possible phase boundary/crossover by 
the horizontal dotted line in Figured]). 

To describe the next three rows in Tabled] note that the order parameters 
(Wo, W4) have a manifestation in gravity. Thus, e.g. if the Zn symmetry 
along t is preserved (Wq = 0), then, due to the 'large A^ volume independence' 
mentioned in Section 12. ![ the free energy F should be independent of /3 



^•^Thc fractional KK scale, related to the 'long string' modes of |2H], can be understood as 
arising from mode shifts of charged fields in the presence of Wilson lines whose eigenvalues 
are uniformly distributed along a circle. 

13 



(equivalently, the action should be hnear in /3). This is satisfied by the 
action in (|6]). Hence for the parameter regime in which the sohtonic D4 
brane is the dominant phase of the theory, the 5 dim SYM dual exists in 
a phase characterized by Wq = 0. Similarly, the non-linear dependence of 
the action (jH]) on L4 implies that the gauge theory dual has W4 7^ 0. These 
results are represented by the 5th row of Table [H One can similarly deduce 
by looking at the (3, L4 dependence of fITOl) that for the phase of the gauge 
theory represented by the black D4 solutioio, W4 = 0, Wq ^ 0, which could 
also be inferred by the t -H- X4 symmetry. This is represented by the 6th row 
of Table [1] In the last row, we have included a phase of SYM5 with (P, AP) 
fermion b.c. along (t, X4) which corresponds to a localized D3 soliton, to be 
discussed in the next section. 

We now come back to the conventional proposal of holographic QCD 
discussed in the beginning of this section; let us define this proposal as a 
combination of three propositions: 

(a) the solitonic D4 brane phase of gravity corresponds to the confine- 
ment phase of YM4, 

(b) the black D4 corresponds to the deconfinement phase of YM4, and 

(c) the SS transition between the two gravity solutions corresponds to the 
confinement /deconfinement transition of YM4. 

In terms of Table [H both the solitonic D4 and the confinement phase of 
YM4 (regarded as a phase of SYM) have Wq = 0, W4 7^ 0; hence both of these 
phases have the same Zn x Zn symmetry and (a) can in principle be valid. 
However, (b) cannot be valid, since the black D4 phase (W^ = 0, Wq 7^ 0) and 



^■^The action (jlOp depends linearly on L4, if we rewrite it by using A5 = A4L4. Hence 
the action per unit length in the X4 direction is independent of L4; therefore, W4 = 0. 
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the deconfinement phase of YM4 (W4 7^ 0, Wq 7^ 0) have different Z^r x Z^ 
symmetry. In particular, the black D4 corresponds to a phase that cannot 
be KK reduced to YM4 at all because of vanishing W4. In fact, it would 
correspond to the Z2 mirror of the confinement phase of YM4 (which can 
be identified with a deconfinement phase of SYM5 which is intrinsically 5 
dimensional). Thus, there must be at least one phase boundary between the 
deconfinement phase of YM4 and the black D4 phase (see Figure [1] for one 
possibility). Regarding (c), if the SS transition persists even at weak cou- 
pling, then it clearly cannot coincide with the confinement /deconfinement 
transition. The reason is that, even at weak coupling the SS transition can 
only occur at /3 = L4 (since even stringy effects should satisfy the Z2 symme- 
try), whereas it is known that the confinement /deconfinement phase transi- 
tion temperature of YM4 goes down (/3c goes up) at smaller values of A4 (see 
Figure [1]) . 

Through the arguments in this section, it is obvious that the black D4 
brane solution does not correspond to the deconfinement phase in the 4 
dimensional Yang-Mills theory and the SS transition does not correspond to 
the confinement /deconfinement transition. In the next section, we propose 
an alternative correspondence, which resolves these problems. 

4 Our proposal for the confinement/deconfinenient 
transition in holographic QCD 

In the previous section, we considered the conventional picture of holographic 
QCD at finite temperature and discussed some of the problems in describing 
phases of YM4. In this section, we propose an alternative picture to resolve 
these problems. 

In the above discussion, we described finite temperature gauge theory by 
choosing AP b.c. for fermions on the brane along the compactified Euclidean 
time. However, let us recall that the gauge theory of interest here is pure Yang 
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Mills theory in four dimensions, which does not have fermions. Fermions 
reappear when the validity condition ([5]) is enforced and Ayjv/ goes above the 
KK scales I/L4. In this sense, fermions are an artifact of the holographic 
method and in the region of validity of the pure YM theory, the periodicity 
of the fermion should not affect the gauge theory results. 

If we start with the 5 dimensional SYM on Si x S}^^ with a periodic tempo- 
ral circle, the partition function is the twisted one: Tr (— l)^exp[— /Ji^sYMs], 
and is not the standard thermal partition function. However, under the limit 
A4 <^ 1 and f3 ^ L4/A4, in which the 4 dimensional YM theory appears, 
the fermions are effectively decoupled; hence F — )■ 0. Thus in this limit the 
'twisted' partition becomes the usual thermal partition function, even with 
the periodic b.c. for the fermions. More precisely, we have 

ZfpYM5) = Tr(-l)^e-'^^sYM5 

= Tre-^^™* + 0{e"'^^^/^^), (A4 ^ 0, L^/X^/S -^ 0). 

7SYM5 _ ry -I3HSYM5 

^(AP,AP) — -"-^^ 

= Tre-'^^Y^" + 0{e-^^^/^*), (A4 ^ 0, L4/A4/3 ^ 0). (13) 

Thus, independent of the boundary condition along the Si, we obtain the 
thermal partition function of YM4. 

As a consequence, it is pertinent to study the phase structure of the 
5 dimensional SYM on 5*1 x S}^^ with (P,AP) boundary condition through 
holographjo. Indeed, in this sector, we will find a phase transition in gravity, 
which continues naturally, at weak coupling, to the confinement /deconfinement 
transition in the 4D Yang-Mills theory. 



^^The phase structure of the Dl branc on a T^ with (P,AP) and (AP,AP) has been 
studied in [27]. Their results are similar to our D4 brane results. The analysis in [29] is 
also related to our study. See also [10] for related comments. 
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4.1 Gregory-Laflamine transition in the gravity theory 
with (P,AP) boundary condition 

In this section, we reconsider the gravity theory described in Section [3l now 
compactified on Si x 5^^ with (P,AP) b.c. for fermions. In this case, contrary 
to the (AP,AP) case, the black D4 brane solution ([H]) is not allowed (see the 
comment below (Q). However a non-trivial phase boundary still exists, as 
we will see shortly. The detailed phase structure for this boundary condition 
is summarized in [T^ ^5\ (with t and X4 exchanged) ; see also [TU] . We will 
now outline the salient features. 

Let us start from large /3. In this region the solitonic D4 brane solution ([2]) 
is thermodynamically dominant. As we decrease (3, it reaches O (^L4/^/\4^, 
which is the order of the effective string length at m = Mq 0, E]- Below 
this value, winding modes of the string wrapping on the /3-cycle would be 
exciteco. Thus the gravity description given by ([2]) would be valid only if 

^ » 4=- (14) 

In order to avoid th is p roblem, we perform a T-duality along the t-cycle and 
go to the IIB framcj, where the solitonic D4 solution becomes a solitonic 
D3 brane uniformly smeared on the dual t-cycle. The metric describing this 



""^^In some previous studies including our works |10[ 136] , the winding modes were not 
correctly evaluated. We would like to thank S. Sugimoto and K. Hashimoto for discussing 
this point. 

^^ In the (AP, AP) case, if we were to perform a T-duality along the i-cycle, the theory 
would be mapped to the OB frame which involves a tachyon. The analysis in the bulk 
gravity in that case would be different from that in the IIB case (we would like to thank 
S. Minwalla for pointing this out to us). However, unlike in the (P, AP) case, the issue 
of light winding modes in the temperature range ([H)) and the consequent need to go to 
a T-dual frame do not arise, since way before we reach temperatures such as P^ . the 
solitonic D4 solution ceases to be the dominant solution, and is replaced by the black D4 
solution at the temperature (jlip . Recall that in the black brane phase the winding number 
is not a conserved quantity, as the time circle is contractible; thus there are no winding 
modes in the spectrum. We hope to come back to a more detailed analysis of the (AP,AP) 
case in a future work. 
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solution is given by 



ds =a 



"\A4A5 



'Tdxl + Uu)dxi\ + ^^ f ^ + dt" + n^dl]^') 



(15) 



Here a't' is the dual of t, hence t' has a periodicity /?' = (27r)^//3 = (27r)^T. 
By considering winding modes on this metric, we can see that the gravity 
description is now valid if /3 <^ //4/A/A4. Thus even when /3 is below (IH1) . 
the gravity analysis in the IIB frame is possible. 

As (3 decreases (hence the dual radius in the IIB frame increases), the 
uniformly smeared solitonic D3 brane becomes thermodynamically unstable 
at a certain temperature. This temperature, called the Gregory-Laflamme 
(GL) instability point, is numerically given by |14p1 : 

i3GL,inst ^ 7.19^. (16) 

It is expected that, even before (3 is lowered all the way to f3GL,inst, the 
smeared solitonic D3 brane becomes meta-stable and undergoes a first or- 
der Gregory-Lafiamme (GL) transition at an inverse temperature Pgl > 
/9GL,msttJ leading to a more stable, solitonic D3 brane localized on the dual 
cycle. This localized solution spontaneously breaks translation symmetry 
along the temporal direction. Although it is difficult to derive the critical 
temperature Pgl precisely, we can approximately evaluate it as follows |29] . 
The metric of the localized solitonic D3 brane is approximately given by that 
of the solitonic D3 on Rg x Sj^^ for a sufficiently large radius /?' of the dual 



^^ Wc obtain this result from Eqn. (12.11) in [H] with p ~ 4 through identification 
i = /3 and f = 1/L4- 

^^Thc GL instabihty was initiahy found in the black string. However similar instability 
exists in the solitonic solutions as well; the reason is easily understood in the Euclidian 
space, where there are no differences between the black objects and solitonic solutions. 



cycle IHEn], 



ds =a' 



u 



VdMP 



t — 1 / 



/3(fi) = !-(?)*. fi„ = V*V^^. .* = ^. (17) 

Here d^ = 2 (see Eq. ([3])). The value of Mq is determined by the smoothness 
condition (as in (jlj)). n is regarded as a radial component of a spherical coor- 
dinate and is related to the cylindrical coordinates (m, t') of ( ITSl) by roughly 
■u^ ~ M^ + t'^. Here we implicitly assume that the soliton is localized at t' = 0. 
This approximation is valid if uq and u are sufficiently small such that we 
can ignore the finite size effect associated with the dual circle [T^ [5U] Pl . 
The classical action of this solution turns out to be 

Ssds/NW, = -C,Lf. (18) 



By comparing this classical action with ([6]) cj, we can estimate that the GL 
transition will happen around 

Pol = ^^ - 8.54^. (19) 

04 A4 A4 

Note that, near this critical point, uq is similar in magnitude to, but less 
than, the dual radius [29], so we expect (fT7|) -( IT9|) to receive some corrections 
from the finite size effect. 

In addition to these two solutions, (fT5|) and (TT7|) . there is another solution: 
solitonic D3 brane non-uniformly smeared on the dual cycle. The metric of 
this solution is perturbatively derived in [13]. This non-uniform solution 
arises at the GL instability point and is always unstable. Although the 



"'^^Later on, in (P5|) . we will consider a situation in which u is large. In that case the 
finite size of the i'-cycle is important and we need to include all the images of the localized 
soliton on this cycle. 

^°Note that the classical action is invariant under the T-duality, hence ([6]) describes 
both the solitonic D4 as well as the smeared solitonic D3. 
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localized SD3 



iion-unifornily smeared SD3 




GL critical point 

1//3CL 



uniformly smeared SD3 
(SD4) 



GL instability point 

^1 PcL.mst 



T=l/0 



Figure 2: The schematic relations among the free energies of various solutions 
in the GL transition (TTOjl [3T| 132] . The red lines denote the stable solutions. 
The blue lines denote the meta-stable solutions. The yellow lines denote 
the unstable solutions. The solitonic D4 (uniformly smeared solitonic D3) 
becomes unstable at the GL instability point f lT6|) . In addition, the localized 
solitonic D3 cannot exist below a certain temperature. 



complete behaviour of the non-uniform solution has not been explored, it is 
expected that it merges with the localized solitonic D3 as shown in Figure [2] 
[I3|[25||3H I32]^H The fact that the GL instability temperature (^E) is higher 
than the approximately obtained GL critical temperature flT^ is consistent 
with this expectation. 

Let us make a few comments on the localized solitonic D3 solution (fT7|) . 
Firstly, through a calculation similar to (|5]), the gravity description can be 
shown to be valid (i.e. stringy effects can be ignored) if A4 ^ P/ L^. Secondly, 
this solution ceases to exist if /3 is too large, as shown in Figure 121 intuitively, 
if /3 were too large, the dual cycle would become smaller than the size of the 
localized soliton, which is not possible. Thirdly, the free energy F = Ssds/P 
of this solution, unlike that of (j2]), is proportional to temperature; hence 
the entropy of this solution is non-zero at 0{N'^) cj and the temporal Zj^ 
symmetry is broken. In addition, because of the non-trivial L4-dependence of 



^^In this paper, we will often use the notation "SDp" for solitonic Dp brane(s). 

^^The entropy derived from the classical action P^ does not depend on temperature. If 
we consider higher order corrections of mq and u in the metric (J17p as in |14| , the corrected 
entropy depends on the temperature. 
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the classical action f llSp we can see that the Zjy symmetry along the X4-cycle is 
also broken. Thus the Polyakov loop Wq and PV4 in this solution are both non- 
zero as shown in Tableland, contrary to the black D4 solution, the localized 
solitonic D3 solution is appropriate for a description of the deconfinement 
phase in the dual gauge theory (which also has Wq 7^ 0,^4 7^ 0). We will 
build on this observation in the next subsection. 

4.2 Gregory-Laflamme transition as a confinement /deconfinement 
transition 

The phase structurcj of the 5 dimensional SYM on Si x S}^^ with the (P,AP) 
boundary condition is shown in Figure |3l The strong coupling region (the 
blue region in Figure [3]) described by type II supergravity and is characterized 
by the GL phase transition which occurs at a temperature given by (TTOl) . In 
the weak coupling region (the green region in Figure [3]) , the 4 dimensional 
Yang-Mills is realized at low temperatures (/3 ^ L4/A4). Although this 
region is common to the (P,AP) phase diagram in Figure |3] and the (AP,AP) 
phase diagram in Figure [H the mirror of this region under the Z2 (/3 ■H- L4) 
does not exist in Figure [3] since the (P,AP) b.c is not Z2-symmetric. 

The main point to emphasize here it that: contrary to the previous phase 
structure in the (AP,AP) case, now the localized solitonic D3 phase has the 
same order parameters (Wq 7^ 0, W4 7^ 0) as the deconfinement phase, thus 
making it plausible that these two phases are smoothly connected. In Figure 
El we have indicated this by assuming the simplest extrapolation through the 
region of the intermediate coupling. 

Our proposal In view of the above observations, we propose a strong 
coupling continuation of weakly coupled 4-dimensional Yang Mills theory as 
shown in Table |2j 



^^In this phase diagram, we do not distinguish the sohtonic D4 and uniformly smeared 
sohtonic D3 because they are essentiahy the same object, ahhough the regions of vahdity 
of their gravity descriptions are different. 
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YM4 



1 ,9 > L4/A4 



confinement 

/ 

'14^0 = 




IIA SUGRA 
(A4 » 1) 
solitonic D4 

(uniformly smeared 
solitonic D3 in IIB) 



Figure 3: The phase structure of five dimensional SYM on Sjj x S}^^ 



with 



the (P,AP) boundary condition. The gravity analysis is valid in the strong 
coupling region (the blue region). The 4 dimensional YM description is valid 
in the green region. The phase structure in the intermediate region, including 
the phase boundary shown by the solid black line, is an extrapolation from 
these two descriptions. 





Weak coupling 


Strong coupling 


Zn X Zn 




4 dimensional YM 


Gravity (P,AP) 




Low temp. 


Confinement phase 


Solitonic D4 phase 


Wo = 0, W4 / 


High temp. 


Deconfinement phase 


Localized solitonic D3 


VFo / 0, W4 / 


Transition 


Confinement/Deconfinement 


Gregory-Laflamme 





Table 2: Our proposal. 
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Of course, our proposal is based on a simple extrapolation between the in- 
termediate coupling regime, and the real story there could be more involved. 
However the mere existence of such a simple extrapolation is a significant 
improvement over the previous proposal in the (AP,AP) case, where we are 
certain that there has to be at least one phase boundary between the de- 
confinement phase and black D4 solution (for the simple reason that their 
order parameters have different values). It is clearly important, therefore, 
to further investigate the nature of the deconfinement phase in terms of the 
localized solitonic D3 solution based on the correspondence outlined above. 

5 New correspondences in holographic QCD 

Using our proposed correspondence, we can explain anew several phenomena 
in the gauge theory from gravity. In this section, we list some of these 
phenomena. 

5.1 Polyakov loop and D3 brane distribution 

In our proposal, we identified the GL transition in the IIB frame to the con- 
finement /deconfinement transition in the gauge theory. In the GL transition, 
the distribution of the D3 branes on the dual circle changes from a uniform 
distribution at low temperatures to a localized one at high temperatures. In 
this subsection, we explain what the corresponding phenomenon is in the 
gauge theory. 

Since the original five-dimensional gauge theory (of the D4 branes) ap- 
pears in the IIA frame, in order to understand the role of the D3 branes, 
we need to consider a T-duality along the temporal direction. Under this T- 
duality, Xq, which is transverse to the D3 brane, is mapped to the gauge po- 
tential Aq on the D4 brane. Thus the D3 brane distribution on the dual circle 
is related to a distribution of ^^'s, which are defined as the eigenvalues of the 



Polyakov loop operator exp(i^fc)5fc; = Pexp(i J A^dt) 



,fc,/ = l,...,Ar. 

kl 
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If the D3 branes are uniformly distributed, the 9^ are also uniformly 
distributed. In that case, by choosing an appropriate gauge, we can take 
9k = 2TTk/N . Then the temporal Polyakov loop operator ([1]) becomes 



1 1 ^ 

— TrPe' Jo' ^odxo = A \^ e'^ik/N ^ q_ ^20) 

AT M ^-^ ^ ' 



N N 

k=l 

Hence this distribution corresponds to the confinement phase. On the other 
hand, if the D3 branes are localized, then 9k are also localized and the 
Polyakov loop becomes non-zero, which characterizes the deconfinement phase. 
This observation is consistent with the entropy arguments in section 14.11 

The above discussion shows a direct relation between the D3 brane distri- 
bution and the eigenvalue distribution of the Polyakov loop operator. This 
latter quantity can sometimes be explicitly evaluated [lOl [331 EH EH [36l [37] . 
E.g. in [in], for four dimensional Yang Mills theory on a small T^, we found 
that the uniform, non- uniform and localized distribution of 9i, all appear, 
with obvious correspondence to similar gravitational solutions. Especially 
the free energies of these three solutions show a "swallow tail" relation simi- 
lar to Figure[2l (See Figure 3 in [ID].) This correspondence strongly supports 
our proposal. 

5.2 Gregory-Laflamme transition as a Hagedorn tran- 
sition 

Our proposal opens up the interesting possibility of a relation between the 
GL transition and the Hagedorn transition. 

It is known that the GL instability is an instability of the KK modes of 
the graviton along the compact circle [I2]. In our case, the KK modes along 
the dual temporal circle, which are associated with the GL instability at (TT6!) 
in the IIB description, are mapped to winding modes around the temporal 
circle through the T-duality [251 [3H]. This then indicates that the type IIB 
GL transition is associated with the excitation of the winding modes of the 
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IIA string^. This phenomenon is similar to the Hagedorn transition in string 
theory [39], where the associated instabihty is caused by temporal winding 
modes. Thus the GL transition in the IIB description might correspond to 
the Hagedorn transition in the IIA description. 

Note that from large N gauge-theoretic calculations also, the usual con- 
finement/deconfinement transition is believed to be related to the Hagedorn 
transition. It has been explicitly shown in weakly coupled gauge theories 
[55| [M]. This makes it plausible that the Hagedorn transition in the Yang- 
Mills theory continues to the Hagedorn transition in the IIA string, which, 
as we argued above, is possibly the dual of the GL transition in the IIB 
supergravity. 

6 Chiral symmetry restoration in Sakai-Sugimoto 
model 

In the previous sections, we have seen that the conventional holographic 
representation of the confinement/deconfinement transition as the SS tran- 
sition is fraught with problems, to circumvent which we proposed in Section 
14.21 a different interpretation in terms of a GL transition. However, the SS 
transition has widely been employed in holographic QCD and purports to 
explain several phenomena in real QCD. In particular, the chiral symmetry 
restoration in the Sakai-Sugimoto model was neatly explained in [11]. In 
their scenario the black D4 brane plays a crucial role. However, since in 
our proposal the black branes do not appear anymore, we need to find an 
alternative idea for realizing chiral symmetry restoration. In this section, we 
discuss how chiral symmetry restoration can happen in the localized solitonic 



^*Note that we have argued in 14.11 that at values of (3 lower than pi)) , winding modes 
get excited. On the other hand, the instabilities from the winding modes arise at the much 
smaller j3 P^ for a large A. It would be interesting to understand the mechanism for the 
separation of these two temperatures in the IIA string theory, although it may be difficult 
since we need to quantize the strings to investigate it. 
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D3 background. 

6.1 Sakai-Sugimoto model and chiral symmetry break- 
ing 

The Sakai-Sugimoto model [H] was proposed to describe low energy liadron 
physics in holographic QCD and elegantly reproduces many aspects of the 
real QCD. We first briefly review this model and show how chiral symmetry 
breaking at low temperatures is realized in terms of the dual gravity. 

The Sakai-Sugimoto model is an extension of the holographic model of 
QCD discussed in section [2]2l Sakai and Sugimoto added, to the A^ D4 brane 
system, Nf D8 and D8 branes which are localized on the x^ circle and fill all 
other directions as follows 

(0) 1 2 3 (4) 5 6 7 8 9 
D^ _____ (21) 

D8/D8 ---- _____ 

Here () denotes the compactified directions. (See Figure H] (a).) This model 
has a free parameter L, which is the asymptotic {u — )■ oo) distance between 
the D8 and D8 brane on the X4 circle. In the original Sakai-Sugimoto model 
[9] L was taken to be L4/2; however, we will let it be general and take values 
in [0, L4/2] (the other half, L e (-^4/2, L4) is related by a reflection and need 
not be considered separately). 

This model has U{Nf) x U{Nf) gauge symmetry on the Nf D8 and D8 
branes, which can be interpreted as a chiral U{Nf)L x U{Nf)fi^ flavor symme- 
try in QCD. As we will see soon, in some situations, the D8 and D8 branes 
merge and the chiral symmetry is broken to a single U{Nf). Sakai and Sug- 
imoto proposed that this is the holographic realization of chiral symmetry 
breaking [9]. 

Let us take a large A^ limit a la Maldacena and, according to the principle 
of holography, replace the A^ D4 branes with a corresponding gravity solu- 
tion. Here we will keep Nf ^ A^ such that we can ignore the back-reaction 
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of the D8/D8 branes onto the background geometry (this is the so-called 
'probe approximation'). In that case, the background (D4 brane) geometry 
is determined thermodynamically (as the dominant classical solution at a 
given temperature, as we have done in sections [3] and Hj) , and the D8/D8 
brane configuration coupled to this background, is determined dynamically, 
with the given distance L as a boundary condition. 

Chiral symmetry breaking in this model happens as follows. At suffi- 
ciently low temperatures, the favoured geometric background is that of soli- 
tonic D4 solution ([2]), as we have seen in section [3] and HI On this background, 
the D8 and D8 cannot extend separately and need to merge as shown in Fig- 
ure H] (b). As a result, the U{Nf) x U{Nf) gauge symmetry on the D8 and 
D8 is broken to U{Nf), representing chiral symmetry breaking (xSB) in the 
dual gauge theory. 

6.2 Chiral symmetry restoration in the black D4 brane 
background 

In the gauge theory, it is expected that chiral symmetry is restored at a 
sufficiently high temperature. Thus if holographic QCD is to work, there 
should be a corresponding phenomena in the dual gravity description. 

In [11], a mechanism for chiral symmetry restoration was suggested by 
considering D8 and D8 in the black D4 brane background ([8]) in the (AP,AP) 
case. Contrary to the solitonic D4 brane geometry (jj]), the u — t plane of the 
black D4 background has a 'cigar' geometry. Thus the D8 and D8 branes can 
wrap the cigar separately as shown in Figure |5] (c). If such a configuration 
is energetically favoured, U{Nf) x U{Nf) gauge symmetry is preserved and 
chiral symmetry is restored. 

In addition to this configuration, another configuration shown in Figure 
|5] (d) is possible in the black D4 background. This configuration is similar 
to Figure |5] (b) and chiral symmetry is broken. In [11], energies of these two 
configurations are compared with the black D4 background. It was found 



27 



(a) D8/D8/D4 set up 



(b) D8/D8 on solitonic D4 background 





XSB 



Figure 4: Brane configuration of the Sakai-Sugimoto model at a finite temper- 
ature. Figure (a) shows the basic brane configuration of the Sakai-Sugimoto 
model before considering the large N limit. The long horizontal direction 



schematically denotes the x ,■ ■ ■ ,x' 

D4 branes (if the D4 branes are localized at x^ 



^ coordinates, which are transverse to the 
■ ■ ■ = x^ = in Carte- 



sian coordinates, the horizontal line can be taken as, for example, x with 
x^ = ■ ■ ■ = x^ = 0). The directions x^,x'^,x^ are common to the D4 and 
D8/D8 branes and are not shown. The directions t and x^ are periodic and 
we should identify 'up' and 'down', and 'front' and 'back' planes in the fig- 
ure, respectively. The square planes on the edges denote T^ = S} x 5*^4 at 
n = oo {x^ = ±oo). The Nf D8 and D8 branes are localized on 5*^4 with the 
asymptotic ( u — t- oo) distance L. The world-volumes of the A^ D4 branes 
are depicted by vertical planes and those of the D8 and D8 by horizontal 
planes. Before taking the large N limit, the D8 and anti-D8 branes extend 
parallelly and chiral symmetry is preserved. In (b) we consider the large N 
limit and replace the N D4 branes by the solitonic D4 brane geometry ([2]). 
This geometry is pinched off at n = uq and D8 and D8 cannot extend to 



the x^, ■ 
broken. 



X- 



direction. Thus they need to merge and chiral symmetry is 
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(c) D8/D8 on black D4 background (d) D8/D8 on black D4 background 




^^^^ D8 



'^ 



D8 




D8 

A 



XSB 



Figure 5: Chiral symmetry restoration in the black D4 brane background ([H]) 
|llj . In this background, although the horizontal direction is still pinched 
off, the t — u plane has a 'cigar' geometry and the D8 and D8 branes can 
separately wrap the cigar as shown in (c). (Recall that the 'front' and the 
'back' of this figure are identical.) In this case, chiral symmetry is restored. 
On the other hand, the configuration shown in (d) is also possible. This 
configuration is similar to (b) and chiral symmetry is broken. 



(see Figure 7 of [TT]) that for L > L^^ = O.I54L4, it is configuration (c) 
which is always favoured, indicating a concurrence of deconfinement and 
chiral symmetry restoration transitions at Tss = '^/Pss = l/-^45 "while for 
L < L^^, a new window Tss < T < 0.154/L opens up where configuration 
(d) is the favoured one, indicating coexistence of deconfinement and broken 
chiral symmetry. 

6.3 Chiral symmetry restoration in the locaUzed soh- 
tonic D3 brane. 

Although the chiral symmetry restoration was explained, as above, in the 
black D4 brane background, our main thesis in this paper is that the black 
D4 solution itself is fraught with problems, if we interpret this solution as the 
deconfinement phase of the four dimensional gauge theory (see Section [3^ . 
The alternative we proposed in Section 14.21 is that the localized solitonic 
D3 brane solution ( flTl) in the (P,AP) b.c. should be taken as the correct 
representation of the deconfinement phase in the dual gauge theory. In view 
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of this, we need to understand chiral symmetry restoration in this solution 
instead of in the black D4 solution. 

One subtlety in the application of our proposal to the Sakai-Sugimoto 
model is the existence of the fundamental quarks from the open string be- 
tween the D4 and D8/D8 brane. These quarks are, of course, not decoupled 
in the 4 dimensional limit A4, L4//3A4 — > 0. Hence, unlike in case of the 
adjoint fermions, where the P and the AP b.c. on the thermal cycle reduce 
to the same quantity in the four dimensional limit (see ( IT3|) ). for the funda- 
mental quarks the P and the AP b.c's differ in the 4D limit, hence, we must 
choose the AP b.c. on the temporal cycle to investigate the thermodynamics 
of the Sakai-Sugimoto model. Since our proposal crucially uses P b.c. for the 
fermions along the temporal cycle, we must address this issue o It turns out 
that with a small addition to our model, viz. that of an imaginary chemical 
potential [ID], we can obtain the desired AP b.c. around the temporal cycle 
for the quarks retaining the P b.c. for the adjoint fermion^. Details of this 
approach will appear in a forthcoming paper [41]. Introduction of such a 
chemical potential does not cause any essential modification to the analysis 
that follows below, as shown in detail in [H] . Consequently, we will proceed 
below with periodic (fundamental as well as adjoint) fermions. 

In our proposal, the gravity analysis has been done in IIB supergravity 
by performing a T-duality along the t-cycle: S^ — ?■ S^, . Thus it is convenient 
to dualize the brane configuration of the Sakai-Sugimoto model to the IIB 
frame. Since all the D branes in the Sakai-Sugimoto model wrap the Si, the 
above T-duality maps the D4 branes and D8/D8 branes to D3 and D7/D7 



^^Of course, for N ^ Nf, the fractional difference in the free energy between the P 
b.c. and the AP b.c, which entirely arises from the fundamental quarks in the 4D limit, 
is 0{Nf/N) — >■ 0. However, such fractional differences are important when we want to 
discuss the free energy of the quarks, which itself is of order 0{Nf/N) compared to that 
of the gluons. 

^^We would like to thank Y. Hidaka for suggesting this point. 
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(e) D7/D7/D3 set up 



(f) D7/D7 on uniform SD3 background 





Figure 6: T-dual of Figure H] along the t-cycle. (e) and (f) are the T-duals 
of (a) and (b) respectively. The distribution of the D3 branes (schematically 
represented by the bunch of vertical lines) on S^, is fixed dynamically. 



branes spreading 

D3_ 

D7/D7 



(0') 1 2 3 (4) 5 6 7 8 9 



(22) 



See Figure [6] (e) also. The distribution along S^, of the D branes in the IIB 
description (related to the gauge field Aq in the IIA description), is deter- 
mined dynamically. In the probe approximation, according to the analyses in 
section H?T| at temperatures below the GL critical temperature f lT9|l . the D3 
branes are distributed uniformly and, above it, they are localized on the S^,. 
The fate of the chiral gauge group in the gravity representation depends on 
the stable configurations of the probe D7/D7 branes in these backgrounds. 

In the uniformly smeared solitonic D3 brane geometry (fT5|) shown in 
Figure |6] (f), the situation is similar to (b) in Figure HI Since the u direction 
is smoothly pinched off at uq, the D7 and D7 have to merge and chiral 
symmetry is broken. 

On the other hand, in the localized solitonic D3 brane geometry, the 
horizontal direction is not fully pinched off. Recall that the geometry (fTTll is 
pinched oS at u = uq, where u ~ y/u"^ + t'"^ and u'^ = {xl + ... +xl)/a''^. Thus 
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(g) D7/D7 on localized SD3 background (h) D7/D7 on localized SD3 background 





xs 



Figure 7: Chiral symmetry restoration in the localized smeared solitonic D3 
brane background. In this geometry, the D3 branes are localized on the S^, 
and the horizontal x^,- ■ ■ ,x^ direction is not fully pinched off. Then we can 
consider two configurations. The chiral symmetry is preserved in (g) and 
broken in (h). 



u can reach zero. As a result, D7 and D7 can extend separately as shown 
in Figure [7] (g) . This configuration is similar to (e) in Figure [6] and would 
restore chiral symmetry. In addition to this configuration, a chiral symmetry 
broken configuration is also possible as shown in Figure [7] (h). 

The last task is the evaluation of the stability of these two configurations 
at temperatures higher than Tql = '^1 Pgl (TT9|) . In the fiat space, the force 
between a single D3 and single D7 (or between a single D3 and a single 
D7) in our configuration is repulsive, since the number of the Neumann- 
Dirichlet open strings between them is 6 |12]. Thus we expect that even 
after taking the large A^ and near horizon limits, the force may be repulsive. 
As a result, the D7/D7 branes, which are separated from the localized SD3 
in the t'-direction, would try to move away as far as possible in this direction. 
However, since the t'-direction is compactified on a circle, the D7/D7 branes 
should then end as being fixed at the point on the t'-circle which is antipodal 
to the localized SD3. E.g. if we put the localized SD3 brane at t' = 0, the 
D7 as well as the D7 branes will both be at t' = /3' /2. See the left diagram 
of Figure [HI In this case, we can effectively restrict the dynamics of these 
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Figure 8: Three classical configurations of the probe D7/D7 brane on the 
localized SD3 background mu — t' and u — X4 space. Green lines correspond 
to the configuration of the D7 and D7 branes in Figure [Tt^g) , which preserve 
chiral symmetry. The red and blue lines correspond to Figure [Tl^h), which 
break chiral symmetry. These three solutions satisfy the same boundary 
condition X4 — )■ L/2 {u — )■ 00). The two x^B solutions appear only for 
L < L^(T). The energy of these solutions are plotted in Figure [9] (a). 



branes to {u — X4) plane (given by t' = (3' /2). Thus the problem reduces 
to finding classical solutions x^i^u) with a boundary condition: X4(m) — )■ L/2 
{u — 7- 00), where we put the D7 at X4 = L/2 and D7 at —L/2. 

The problem stated above is difficult to solve precisely near the GL tran- 
sition since the background metric ( !T7|) is only an approximate description. 
The metric around the D7/D7 brane (at t' = (3'/2) becomes more and more 
accurate, however, when T ^ Tgl- As a result, we solve for the stable con- 
figuration in this limit. Details of the calculation are presented in appendix 
Rl We find three solutions, corresponding to the D brane configurations of 
Figures [7] (g) and (h), as shown in Figure [HI In the appendix, we compare 
the classical DBI actions of these solutions numerically, as depicted in Figure 
|9](a) (see also the phase diagram in Figure [TOj) : as one increases temperature 
beyond a certain value, chiral symmetry is restored, the transition being of 
first order. Therefore, one can see that even in our proposal, similarly to the 
case of the black D4 [H] , we can explain chiral symmetry restoration at high 
enough temperatures. 
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(a) Sdbi vs. L 



(b) L vs. wi 



Figure 9: (a) The DBI action for probe D7 branes vs. L. We evaluated ( 125|) 
for the three solutions (one X'S' and two X'S'-B solutions). Since these DBI 
actions diverge, we regularize them by subtracting the DBI action for the xS 
solution. A first order transition occurs at L = L^. (b) The relation between 
Ui and L in the xSB solution f l27|) . Two x^B solutions exist for a given L 
ii L < L^. The colour coding in these figures is the same as that in Figure El 
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deconfinement x.S'i? 



Tgl 
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Figure 10: Preliminary phase diagram of the Sakai-Sugimoto model. We 
numerically evaluate Lc{T) by using the high temperature approximation 
and the result near Tgl is not reliable. The T coordinate is scaled by A4/L4. 
In higher temperature, the xS phase is favoured. 
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7 Conclusions 

In this paper, we showed that the conventional representation of the con- 
finement /deconfinement transition in holographic QCD has several problems 
and proposed an alternative representation which resolves these problems. 

As mentioned earlier, problems similar to the above had also been en- 
countered in the study of two dimensional bosonic gauge theory in [TU] . 
This indicates that the issues addressed in this paper are rather general 
in the discussion of holography for non-supersymmetric gauge theories at fi- 
nite temperatures. To elaborate, in the standard holographic procedure, a 
p-dimensional non-supersymmetric gauge theory is first constructed through 
the KK reduction of a {p + l)-dimensional super Yang-Mills theory on a 
Scherk-Schwarz circle. The {p + l)-dimensional SYM at large A^ can be 
mapped to a scaling limit of Dp brane geometries [2] . At finite temperatures, 
because of the two compact cycles (temporal and Scherk-Schwarz), several 
distinct solutions (depending on boundary conditions) appear in gravity as 
shown in section |3] and HI solitonic Dp (equivalently, uniformly smeared soli- 
tonic D(p — 1)), localized solitonic D(p — 1), and black Dp. The black Dp 
brane solution appears at high temperatures in the (AP,AP) case, while the 
localized solitonic D{p — 1) brane solution is the high temperature phase 
in the (P,AP) case. A table similar to Table [U can again be constructed, 
where the appropriate order parameters would be Wp, Wq; these would again 
appear to favour the localized solitonic D{p — 1) phase as the more suitable 
representation of the deconfinement of p- dimensional YM theory (rather than 
the more conventional black Dp phase which appears only in the (AP,AP) 
b.c). Following this logic, a p-dimensional analogue of our proposal (see 
Section 14. 2p would appear to give a better description of holographic QCD 
in p dimensions. In particular, we believe that the Gregory-Lafiamme tran- 
sition between the solitonic Dp and localized solitonic D{p — 1), would, as 
in this paper, be related to the confinement /deconfinement transition in the 
p-dimensional gauge theory. 
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7.1 Further questions 

In order to further understand holographic QCD through the above proposal, 
it would be of interest to address the following questions. 

Transition temperature From (IT9]) . the critical temperature of the con- 
finement /deconfinement transition (the GL transition) would be 0{\4,/Li). 
On the other hand, holographic QCD seems to predict that the square root 
of the QCD string tension is 0{y/\l/ L^), as in (IT^ . whereas the glueball 
masses are 0(1/L4) [6]. It would be important to understand the reason for 
the separation of these scales and how they evolve from strong coupling to 
weak coupling. 

Quantitative correspondence Although our new proposal reproduces 
the known qualitative features of 4 dimensional Yang-Mills theory, it does 
not automatically lead to a quantitative agreement. For example, the free 
energy of the deconfinement phase of the YM theory at a sufficiently high 
temperature must be proportional to T^, since the coupling becomes weak 
and the theory becomes approximately conformal. However the free energy 
of the localized solitonic D3 is proportional to T (as can be seen from ( IT8|) ). 
This is not entirely surprising since the functional form of F{T) can change 
as one evolves from weak coupling to strong coupling. Furthermore, one 
has to exercise caution in defining a high temperature limit of YM4 in the 
holographic context since the temperature must always remain much smaller 
than the KK scale. 

YM4 from SYM5 with (AP,AP) In this article, we have emphasized 
the correspondence between YM4 and SYM5 with (P,AP) b.c. However, as 
we pointed out in flT3|) . the 5 dimensional SYM with (AP,AP) b.c. should 
also be related to YM4, since the boundary condition becomes irrelevant in 
the limit A4 ^ 1 and /3 ^ L4/A4. A possible way this correspondence may 
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work is as follows. In the (AP,AP) case, if we treat the black D4 solution 
(the blue region in the bottom half of Figure [1]), which is not related to YM4, 
as irrelevant, and focus on the solitonic D4 brane, its winding modes would 
appear to be light around the temperature (1141) . A Hagedorn transition 
parallel to the (P,AP) case might occur above this temperature and might 
continue to the confinement /deconfinement transition in YM4. However, in 
order to investigate it further through gravity, we need to understand the 
gauge/gravity correspondence in the OB frame as mentioned in footnote dSi 

Real time It would be important to explore what geometry corresponds to 
the deconfinement phase in the real time formalism. An understanding of this 
would allow us to address dynamical properties of the deconfinement phase, 
e.g. transport properties. Since previous results in this area were based 
on the black brane solutions, it would be important to see how well-known 
results such as the viscosity bound [13] can be derived in our proposao. 

AdS/CMT and chemical potential dependence The SS transition for 
D3 branes has been studied in the context of the AdS/CMT correspondence 
to investigate the superconductor/insulator transition in 2+1 dimension [il| 
H5] . In these studies, a chemical potential for a U{1) charge was introduced 
and the phase structure involving this chemical potential has been derived. 
It would be interesting to ask whether our proposal has any bearing on these 
studies, e.g., whether the GL transition analogous to the one discussed here 
can be a candidate for the superconductor /insulator or some other transition. 
A possible line of investigation could be to study the chemical potential 
dependence of Tc for the GL transition vis-a-vis that for the SS transition 
and see whether any qualitative differences appeac^. 



^^We would like to thank E. Kiritsis for emphasizing this issue to us. 
^®The chemical potential dependence for the confinement /deconfinement transition di- 
rectly from gauge theory has been studied in a one dimensional model j37j . 



37 



Acknowledgement 

We would like to thank Avinash Dhar, Saumen Datta, Koji Hashimoto, 
Yoshimasa Hidaka, Elias Kiritsis, Matthew Lippert, Shiraz Minwalla, Vasilis 
Niarchos, Tadakatsu Sakai, Shigeki Sugimoto and Tadashi Takayanagi. We 
would also like to thank Ofer Aharony, Matthew Lippert, Shiraz Minwalla, 
Rob Myers and Vasilis Niarchos for a critical reading of the manuscript and 
for offering suggestions for improvement. We would like to thank the organiz- 
ers of the sixth regional String theory meeting in Milos, Greece (June 2011) 
for hospitality during the final stages of this work and for an opportunity to 
present this work. T.M. would like to thank KEK for hospitality where part 
of the work was done. T.M. is partially supported by Regional Potential 
program of the E.U. FP7-REGPOT-2008-1: CreteHEPCosmo-228644 and 
by Marie Curie contract PIRG06-GA-2009-256487. 

A Probe D7 brane solution in T ^ Tgl 



In this appendix, we show the calculation of the probe D7/D7 solution in the 
localized solitonic D3 brane background (fTTl) in the T ^ Tgl limit. 



As we argued in section 16. 3[ the D7/D7 would be localized at t' = /3'/2. 
Then, if the temperature is sufficiently high, we can use the Newtonian ap- 
proximation [131 [301 SHI SZ] around the position of the D7/D7 and the back- 
ground metric in this region would be approximately described by 

3 



ds' 
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Here the sums in F and $ are contributions of the mirror images of the 
locahzed SD3 on the t'-cycle and we have assumed $ ^ 1. We consider the 
classical solution of the the D7/D7 brane on this background. Since these 
branes would be localized at t' = /3'/2, the non-trivial configuration of the 
branes would be X4(m) only. Then the induced metric on the D7/D7 becomes 

3 



(isf^y = a' 



+ (Fo~'/' (1 + 2<l>o) x^iuf + fJ/') du' 

_ XiLiP'^ sinhx-x TT^X^ Pu 

^° - ^'*'=^'/2 - 64vr4 x3cosh^(x/2)' ^° " ^2^^^°' "" " 2^" ^^^> 



Here we have used ds = 2 and (3 = (27r)^//3'. Now the DBI actioio of the 
D7/D7 is given as 

Sdbi =NfTj I d^xe-'^^ydetgD7 oc /" rfn (1 - $0) u^^/il + 2%) x^{u)^ + Fq. 

(25) 
and the equation of motion for X4{u) from this action is 

^ / {l + %)u'x^{u) \ ^ ^ 
du \^{l + 2%)x^{uf + Fo) 

Recall that the boundary condition for Xi{u) is X4 = L/2 {u — )■ 00). (From 
now we restrict X4^{u) > by assuming that the negative solution is obtained 
by -X4{u).) 

From this equation, we can find two classical solutions. One is a constant 
solution: X4{u) = L/2. This solution corresponds to the xS line in Figures 
[H]and[n] (a). The other solution is described by 



x,{u)= I du'il-Mu))\l ^^,^;^y_^ . (27) 



Ml 



"^^ We can verify the localization of the D7/D7 brane at t' = /3'/2 by evaluating the 
fluctuation of t' in this action. 
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Here we have assumed the existence of a turning point Mi, which satisfies 
X4,{ui) = and X4(ui) = oo. This solution corresponds to the xSB profile 
in Figures |8] and [91 The value of mi can be determined by the condition 
X4{u = oo) = L/2. We numericall}o find the relation between ui and L, 
which is plotted in Figure M (b). There is a critical value L^:{T), and, if L 
exceeds it, no x^B solution exists whereas two solutions exist if L < L*. 
One solution (the blue lines in Figure [8] and |9]) is always unstable compared 
to another one (the red lines). 

By comparing the DBI action ( 125|) for these two x^B solutions and the 
xS solution, we find that the xS solution is stable for larger L. A first order 
transition occurs at L = Lc{T) and the x^B solution becomes stable below it. 
See Figures[n](a) andfTUl Note that Lc{T) and L^:(T) are decreasing functions 
of temperature. Thus the xS phase is favoured at higher temperature, which 
is consistent with QCD. 

One important question is whether chiral symmetry is restored at T = 
Tgl- In our numerical calculation, we find that, if L < Lc{Tgl) = O.I2L4, the 
confinement/ deconfinement transition and the chiral symmetry restoration 
transition are separated as depicted in Figure fTOl However, we have used the 
high temperature approximation in this calculation and the result would not 
be reliable near Tgl- In addition, the background metric (fTTj) was derived 
using an approximation which also does not work near this temperature. In 
order to improve the result near Tgl, we may need precise numerical study 
as in l3n. 
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